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ABSTRACT 

We suggest that the intrinsic, stellar initial mass function (IMF) follows a power-law slope 7 = 2, 
inherited from hierarchical fragmentation of molecular clouds into clumps and clumps into stars. 
The well-known, logarithmic Salpeter slope T = 1.35 in clusters is then the aggregate slope for all 
the star-forming clumps contributing to an individual cluster, and it is steeper than the intrinsic 
slope within individual clumps because the smallest star-forming clumps contributing to any given 
cluster are unable to form the highest-mass stars. Our Monte Carlo simulations demonstrate that 
the Salpeter power-law index is the limiting value obtained for the cluster IMF when the lower-mass 
limits for allowed stellar masses and star- forming clumps are effectively equal, m\ = M\ a . This 
condition indeed is imposed for the high-mass IMF tail by the turn-over at the characteristic value 
m r . ~ 1 M Q . IMF slopes of V ~ 2 are obtained if the stellar and clump upper-mass limits are also equal 



'up 



= M up ~ 100 M Q , and so our model explains the observed range of IMF slopes between r ~ 1 to 
2. Flatter slopes of T — 1 are expected when M\ > m up , which is a plausible condition in starbursts, 
where such slopes are suggested to occur. While this model is a simplistic parameterization of the 
star-formation process, it seems likely to capture the essential elements that generate the Salpeter tail 
of the IMF for massive stars. These principles also likely explain the IGIMF effect seen in low-density 
star-forming environments. 

Subject headings: stars: formation — stars: luminosity function, mass function — stars: massive — 
stars: statistics — galaxies: star clusters: general — galaxies: stellar content 



1. INTRODUCTION 

The mass distribution of stars at birth, known as the 
initial mass function (IMF), is perhaps the most impor- 
tant fundamental parameterization of the star formation 
process (e.g., Lada 2009; McKee & Ostriker 2007; Adams 
2003). When Salpeter (1955) first empirically evaluated 
the IMF, he identified a simple power-law distribution 
for stellar mass m, having exponent 7 = 2.35: 



n(m) dm cx m 7 dm 



(1) 



While it is now clear that the IMF has a characteristic 
mass below which its form turns over, yielding additional 
distinct regimes at lower masses (e.g., Kroupa 2001; 
Chabrier 2003), Salpeter's original prescription never- 
theless has remained surprisingly robust for m > 1 M Q . 
Since the stellar light from galaxies is dominated by these 
more luminous stars, the "Salpeter IMF" has proven vi- 
tal to the entire field of galaxy evolution, in addition to 
star formation and stellar populations. 

Yet, more than half a century after Salpeter's (1955) 
seminal work, the physical factors yielding the value of 
7 = 2.35 for the power-law index remain elusive and 
poorly understood. The origin of the IMF continues to 
be a topic of intense discussion, and the reader is re- 
ferred to recent reviews by Lada (2009), Clarke (2009), 
Elmegreen (2009), McKee & Ostriker (2007), Bonnell, 
Larson & Zinnecker (2007) for comprehensive overviews 
on IMF theory. 

One obstacle to gaining deeper physical understanding 
is the fact that the Salpeter slope emerges from a vari- 
ety of simulations dominated by different mechanisms. 
In particular, both the core collapse (e.g., Krumholz 
et al. 2010) and competitive accretion scenarios (e.g., 



Bonnell et al. 2003) are able to reproduce the Salpeter 
slope. Effects that are explored in modern simulations 
include decaying vs driven turbulence, isothermal vs non- 
isothermal equations of state, and inclusion or exclusion 
of magnetic fields, among other factors; the IMF slope 
generally does not provide a strong discriminant on this 
large parameter space (e.g., Clarke 2009; McKee & Os- 
triker 2007). 

Here, we suggest an explanation of the Salpeter slope 
that may provide the simplest basic physical understand- 
ing, but which may correspond to a framework within 
which the more specific physics of actual star formation 
operates. 

2. FROM THE CLUSTER MASS FUNCTION TO THE IMF 

Most, if not all, stars form in clusters (e.g., McKee & 
Ostriker 2007; Lada & Lada 2003), which in turn form 
with an initial cluster mass function (ICMF). Like the 
stellar IMF, the cluster ICMF is also described well by a 
power-law distribution in cluster mass M: 



N(M) dM oc M - ' 9 dM 



(2) 



Again like the stellar IMF, the ICMF power law index is 
also found to be robust and largely invariant, but slightly 
shallower, fj ~ 2.0, than the Salpeter slope, based on 
clusters covering a wide range of scales, including super 
star clusters (e.g., Meurer et al. 1995), massive and open 
clusters (e.g., Elmegreen & Efrcmov 1997; Zhang & Fall 
1999; de Grijs et al. 2003), sparse OB groups (Oey et al. 
2004), and the Hn region luminosity function (e.g., Oey 

6 Clarke 1998; Kennicutt et al. 1989). 

The similarity in the slopes of the IMF and ICMF, 

7 = 2.35 vs /3 = 2.0, respectively, has been previously 
noted (e.g., Elmegreen 2006). Since star formation is a 
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hierarchical process, with smaller units fragmenting from 
larger ones, we can examine the origin of the IMF as 
the relationship between the IMF and ICMF. The ICMF 
slope (5 = 2 is a value that seems simpler to understand. 
As is often pointed out, this power-law index corresponds 
to a uniform distribution of power between large and 
small scales (e.g., Elmegreen 2006). Furthermore, Zin- 
necker (1982) showed that an IMF slope 7 = 2 results 
simply from Bondi-Hoyle accretion m oc m 2 among the 
protostellar masses, with the mass function n(m, t) evolv- 
ing from some arbitrary initial distribution: 



dn 
~dt 



d ,. . 
+ —{mn) = 
dm 



(3) 



This yields a simple power law IMF having 7 = 2 (see 
also Bonnell et al. 2007; Hsu et al. 2010). Thus 7 = 2 
is a value that is intuitive, with simple and reasonable 
physical bases from both the parent ICMF equipartition 
arguments and build-up from simple accretion. So then 
why is the Salpeter slope of the IMF slightly steeper than 
7 = 2? 

Elmegreen (2009) summarizes the three approaches to 
generating the IMF: 1) Fragmentation, in which the core 
mass function and IMF is a direct mapping from the par- 
ent cloud fragmentation; 2) Accretion, in which the IMF 
is the product of protostar accretion processes not neces- 
sarily linked to how clouds fragment; and 3) Interruption 
of 1) or 2), in which the form of the IMF is also deter- 
mined by factors that preferentially limit fragmentation 
or accretion for high or low-mass stars. If we adopt our 
argument above that, simplistically, the IMF should in- 
herit 7 = 2 from the ICMF according to both 1) and 2), 
then 3) is a likely candidate to explain the deviation of 
the Salpeter slope: there is likely some process that lim- 
its the formation of high-mass stars in favor of low-mass 
ones. 

In a fully hierarchical star formation scenario, the 
masses of the parent clouds, and hence, cluster masses, 
are determined before the masses of the star-forming 
clumps, and in turn, the constituent stars. If the small- 
est clumps have masses on the order of stellar masses, 
then the formation of high-mass stars will necessarily be 
limited in such clumps, relative to the low-mass stars. 
In recent years, it has been proposed that the integrated 
galaxy IMF (IGIMF) has a steeper slope than the IMF 
because the most massive stars may not form in smallest 
clusters (Kroupa & Weidner 2003; Weidner & Kroupa 
2005, 2006; see also Elmegreen 1999), causing the max- 
imum stellar mass m up to depend on the cluster mass. 
This would cause the IGIMF to favor lower-mass stars. 

We propose that the difference between the Salpeter 
slope 7 = 2.35 and the ICMF slope (3 = 2 is caused by 
a similar situation governing the distribution of stellar 
masses within individual clusters. If, following a fully 
hierarchical scenario, stellar subgroups are formed out 
of clumps with pre-detcrmined masses, and if the clump 
mass distribution also generates a dependence of m up on 
clump mass M, then the resulting IMF for the entire clus- 
ter will be slightly steeper than the clump mass function. 
This scenario has essentially been proposed by Elmegreen 
(1997), who presented this as random sampling from a 
fractal mass distribution. Here, we present a parame- 
terization that is much simpler than Elmegreen's model, 



and we show that the Salpeter slope simply results from 
the overlap between the stellar mass range and the clump 
mass range. 

3. RELATIVE MASS RANGES OF CLUMPS AND STARS 

We model the generation of stars in a cluster as the 
sum of all stars forming from individual clumps in the 
parent molecular cloud. We assume that the hierarchi- 
cal fragmentation of the cloud into clumps follows the 
same (3 — 2 power-law distribution found for the ICMF, 
and further, that the fragmentation of clumps into stars 
does the same, arguing that self-similar physical mecha- 
nisms govern all these fragmentation processes on smaller 
scales. This breaks down around the characteristic value 
of m c ~ 1 M Q , suggested to be linked to the Jeans mass 
(e.g., McKee & Ostriker 2007; Clarke 2009). Since we 
are interested only in the upper IMF, however, scale-free 
cloud fragmentation plausibly dominates in this regime. 
But we note that our model contrasts with those in which 
the IMF slope results from physics that causes fragmen- 
tation into a steeper core mass distribution (e.g., Padoan 
& Nordlund 2002). 

We construct Monte Carlo simulations drawing the 
clump mass function from the j3 = 2 power law, and 
then the stellar mass function within each clump from 
the same distribution, within lower and upper stellar 
mass limits m\ = and m up . The entire cluster is then 
the aggregate of all the stars formed in all the clumps, 
and the IMF is the composite for this aggregate. Our 
default sampling algorithm allows rejection of stars that 
cause the clump mass limit to be exceeded, and contin- 
ues sampling until the total stellar mass is within m\ Q of 
the specified clump mass. This is similar to, but a bit 
more strongly sorted, than the "sorted sampling" algo- 
rithm of Weidner & Kroupa (2006; see also Elmegreen 
2006; Parker & Goodwin 2007). 
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Figure 1. Dependence of fitted logarithmic IMF slope T on the al- 
lowed mass range for stars {m\ Q to m up ) relative to that for clumps 
{M\ Q to M up ). Panels a - c are based on sorted sampling, and 
panel d is based on unsorted sampling. To ease comparison, the 
black symbol shows the same model parameters, and the dotted 
line shows the Salpeter slope value. 

Figure [T] shows how the observed logarithmic IMF 
slope r, determined from fits weighted by the inverse 



Origin of Salpeter Slope 



3 



of the poisson errors, depends on the mass ranges for 
the stars relative to the clumps. If the clump masses 
are all larger than the allowed stellar masses, then the 
stellar masses are unconstrained, and the parent (3 — 2 
distribution is reproduced, translating to T = 1 on a 
logarithmic plot. This is apparent in Figure QJi, hav- 
ing a fixed allowed stellar mass range of 1 to 100 M , 
and clump upper-mass limit M up = 10 4 M Q . Figure [TJi 
shows r as a function of the lower clump limit M\ , and 
as M\ — > m up , r — > 1. But as allowed cloud masses de- 
crease to values below m up , the IMF slope steepens, since 
formation of the most massive stars is excluded. We 
see that as M\ a — > mi OJ the value approaches T = 1.4, 
essentially the Salpeter value (black point). Figure [TJ) 
confirms that this effect is driven by the allowed clump 
lower-mass limit, rather than the upper limits on either 
the clumps or stars. This model allows both M up and 
m up to be essentially limitless, 10 6 M Q ; the results are 
almost identical to Figure [TJz in approaching the same 
Salpeter slope as M\ a — > m\ Q . This steepening effect be- 
tween [3 — 2 and the Salpeter value can be seen in models 
by Elmegreen (1997) and Weidner & Kroupa (2006). In- 
deed, Elmegreen's model allowed an effective equivalent 
M\ ~ 1 Mq, and Weidner & Kroupa adopted an analo- 
gous M\ Q = 5 M Q ; both of these clump lower- mass limits 
are well within the regime where the IMF steepening is 
expected, according to our simulations. 

IMF slopes steeper than the Salpeter value can be ob- 
tained if both M\ Q <C m up and M up approaches stellar 
values. Figure [TJ: assumes M\ a = mio, and shows T as 
a function of M up . For M up ^ m up , T ~ 1.4, as found 
before for the assumed M\ Q — m\ . As the clump upper- 
mass limit is decreased, T steepens, reaching T = 1.9 
for M up = m up = 100 Mq. The results are essentially 
identical even if m up is limitless, although the maximum 
stellar mass will be limited by M np at low M. Weidner & 
Kroupa (2006) also show that the degree to which the ag- 
gregate IMF steepens depends slightly on the sampling 
algorithm used. A "sorted sampling" algorithm, as we 
use, allows continued sampling of lower-mass stars after 
the highest-mass stars can no longer fit within the alloted 
mass, and naturally induces slightly steeper slopes than 
algorithms that stop sampling as soon as any drawn star 
causes the clump mass limit to be reached. Figure [TJf is 
the same as panel c, but with unsorted sampling, con- 
firming that the sorting scheme is not a strong effect. 

Thus we see that the slope steepening beyond the input 
(3 = 2 value is driven by allowing clump masses <C m up . 
In these lowest-mass clumps, the highest-mass stars can- 
not form, thus causing the aggregate cluster to slightly fa- 
vor lower-mass stars beyond the conventional IMF. From 
a physical standpoint, this simply means that clumps too 
small to produce the highest-mass stars can still produce 
low-mass stars. Furthermore, the Salpeter IMF results 
when M\ Q ~ m\ while randomly sampling both stars and 
parent clumps from the same (3 — 2 power-law distribu- 
tion. Thus ironically, the critical parameter for steep- 
ening the IMF and the IGIMF is not M up , but rather 
M lo . 

Elmegreen (1997, 1999) developed a model for IMF 
generation based on random sampling of stellar masses 
from hierarchical fractal structure. Our model distills 
the effects seen in Elmegreen's model to the most sim- 



plististic level, and it more directly captures the domi- 
nant effect seen in that work, revealing M\ as the critical 
parameter. 

4. DISCUSSION 

The essential condition that drives the steepening of 
the IMF from 7 = 2 is that M\ D <C m up . This need not 
be interpreted literally, but the smallest clumps must be 
incapable of producing the highest-mass stars; our model 
is independent of star-formation efficiency, provided that 
it is essentially constant. To attain the Salpeter value, 
in particular, requires the equivalent of M\ Q ~ m\ . Be- 
cause the full IMF flattens strongly near a characteris- 
tic value m c ~ 1 M©, this limits the power-law form of 
the high-mass IMF to apply only above this characteris- 
tic mass, thus m\ ~ m c , recalling that here we consider 
only the upper, Salpeter IMF. Furthermore, because stel- 
lar objects continue to form with masses 1-2 orders of 
magnitude lower, clump fragmentation clearly continues 
to masses below m c as well, and so we know that the 
relevant lowest-mass clumps having m\ ~ m c do form 
stars. Thus, for the power-law upper IMF, M\ Q < m\ Q , 
which is exactly the condition needed to induce steepen- 
ing of the aggregate cluster IMF to the Salpeter value 
in the power-law regime above m c . In other words, m c 
acts as an effective lower cutoff to the IMF power-law 
distribution for both clump and stellar masses, render- 
ing Mio = mi . Indeed, C ls O observations of M17-SW 
show that molecular clump masses follow a power law 
of T = -0.72 ± 0.15 down to at least 1 M Q (Stutzki k 
Glisten 1990), while the observed IMF in the closely as- 
sociated region M17-SWex shows T = 1.3 ± 0.2 (Povich 
& Whitney 2010). 

If clumps fragment according to a (3 — 2 power law, 
then the range of observed IMF slopes should generally 
not fall below this value. Similarly, steeper IMF slopes 
reach T - 2 for M up = m np = 100 M , when M Xo ~ 
mio- A number of authors have shown compilations of 
measured cluster slopes, and 1 < T < 2 does appear to 
be a fairly well-defined range of allowed slopes for the 
high- mass IMF (e.g., Elmegreen 1999, Fig. 1; Bastian 
et al. 2010, Fig. 2): what appears to be a large scatter 
about r = 1.35 (Elmegreen 1999) can also be interpreted 
as variation within the range allowed by our model. 

The basic premise of our model is that clumps frag- 
ment from clouds according to the (3 = 2 power law, and 
in particular, that stars formed within each individual 
clump are still formed with a 7 = 2 power law. Thus 
if newly-formed stars can be identified with their na- 
tal sub-groups in the youngest star-forming regions, the 
IMF should appear to be closer to T = 1 in such groups. 
It would be interesting to evaluate the IMF slopes for 
sub-groups within massive star-forming regions, and to 
compare these to the aggregate cluster IMF. Also, obser- 
vations focusing on specific subregions within extremely 
young clusters might be more likely to obtain IMF slopes 
closer to T = 1, while more complete measurements for 
entire clusters might obtain the Salpeter value or higher. 
These are difficult tests, since dynamical evolution will 
quickly mix stars born from different clumps. 

It is often suggested (e.g., Elmegreen 2004, 2009) 
that the IMF is slightly flatter in starbursts and rich, 
young, super star clusters. In our model, recovering 
an IMF slope of 7 = 2 occurs naturally if the clump 
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mass M\ > m up . The Arches cluster near the Galactic 
Center is the best-known example of a flat IMF, with 
the most recent measurements of the high-mass IMF of 
r = 0.91 ± 0.08 (Kim et al. 2006) and T = 1.1 ± 0.2 (Es- 
pinoza et al. 2009); these values are consistent with our 
model. Regions of extreme star formation are expected 
to have higher thermal Jeans masses owing to radiative 
feedback from massive stars (e.g., Larson 2005; Murray 
2009) therefore raising the minimum clump mass M\ a . 
Our slope flattening would only occur if the clump mass 
range is increased without a commensurate increase in 
the allowed stellar mass range; and indeed, m up appears 
to be remarkably constant and independent of star for- 
mation environment (e.g., Oey & Clarke 2005; Weidncr 
& Kroupa 2004). 

If the smallest clusters form out of single clumps, with 
no further fragmentation, then the cluster IMF is the 
same as the clump IMF and there is no steepening of 
the aggregate slope. Hence, such single-clump clusters 
should show the parent 7 = 2 stellar IMF. Further, if the 
smallest clusters do approach limiting masses that are on 
the order of the largest clump masses, this should cause 
an analogous steepening of the cluster mass function un- 
der these circumstances. As described above, the ICMF 
generally has a slope /3 = 2, but it would be interesting 
to evaluate any trends. We also note that our model docs 
not necessarily imply an IGIMF steepening itself: smaller 
clusters can still fragment into clumps yielding the full 
stellar mass range for all, including the smallest, clumps. 
So if the smallest clusters are still capable of forming 
stars up to masses m up , as suggested by the results of 
Lamb et al. (2010), then no IGIMF steepening will oc- 
cur. Thus again, the lower clump mass range relative 
to the upper stellar mass range determines the steepen- 
ing for both the cluster stellar IMF and the IGIMF. We 
do however suggest that the IGIMF effect seen in low- 
density environments (e.g., Lee et al. 2009; Hoversten & 
Glazebrook 2008) is real and occurs because the smallest 
clusters have masses below m up . 

We stress that our model is a simplistic parameter- 
ization of hierarchical star formation, but it may cap- 
ture essential elements that govern the generation of the 
Salpeter tail of the IMF. In particular, our model is based 
on a scenario that is fully hierarchical, with clouds frag- 
menting into clumps and clumps into stars, with each 
fragmentation process based on a f3 = 7 = 2 power-law 
mass distribution. "Fragmentation" need not be literal, 
but simply means that the ultimate mass apportionment 
from the parent to descendant units follows this power- 
law distribution, whether by competitive accretion or 
by well-defined fragmentation. The steepening of the 
IMF to the Salpeter value is then induced by the inabil- 
ity of the lowest-mass clumps to form the highest-mass 
stars. Maschberger et al. (2010) have carried out a star- 
formation simulation of a large volume, which follows 
the hierarchical fragmentation of ISM into clusters and 
down to core-like sink particles. Their results are consis- 
tent with our scenario, showing that clusters are mergers 
of products from multiple clumps, and that the aggre- 
gate IMF in such mergers is slightly steeper than in the 
individual merging units. 



5. CONCLUSION 

We propose that the empirical, logarithmic T = 1.35 
Salpeter power-law slope of the IMF for high-mass stars 
originates from a universal, linear 7 = 2 IMF slope for 
stellar sub-groups formed from clumps within individual 
clusters. Our Monte Carlo simulations show that T = 1.4 
is the limiting value of the IMF slope when the lower- 
mass limit for both clumps and stars is equal, M\ = m\ -, 
this condition holds for the high-mass tail of the observed 
IMF, since it turns over below the characteristic mass 
m c ~ 1 M Q , effectively setting the lower-mass limit to 
both TOio and M\ . The steepening to the Salpeter value 
occurs because the highest-mass stars cannot form in the 
lowest-mass clumps. Thus, a critical factor to examine 
in star formation theory is the effective mass range for 
clumps relative to the output stellar masses. Our model 
is analogous to the mechanism for steepening the IGIMF 
proposed by Kroupa & Weidner (2003), but a universal 
Salpeter slope does not necessarily imply that the IGIMF 
should steepen to even higher values, since the smallest 
clusters still may be capable of forming the highest-mass 
stars. However, we stress that the same principles involv- 
ing the relative mass ranges of the parent and descendant 
units also applies to the IGIMF and may be responsible 
for the apparent steepening in regimes dominated by the 
smallest-scale star formation (e.g., Lee et al. 2009). 

Our model is based on a fully hierarchical scenario in 
which the physics of fragmentation from large scales to 
individual stars takes place in a self-similar manner, all 
with a -2 power-law mass distribution. Thus, this does 
not require literal fragmentation, but simply that the fi- 
nal mass distribution among both the parent and descen- 
dant units is distributed accordingly. With this simple 
condition, we show that the aggregate high-mass IMF 
slope is limited between values of T = 1 and 2, which is 
consistent with observations. Flatter values near T = 1 
occur when M\ > m up , which may explain the suggested 
IMF flattening for regions of extremely intense, massive 
star formation. While our model is only a rough param- 
eterization of the star formation process, we suggest that 
it may capture the fundamental effects that generate the 
Salpeter slope of IMF for high-mass stars. 

Thanks to Cathie Clarke, Joel Lamb, and Fred Adams 
for useful discussions and comments, and to the referee, 
Michael Reid. I gratefully acknowledge the Institute of 
Astronomy, Cambridge for visitor support, and funding 
by NSF grant AST-0907758. 
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